Bond-operators and triplon analysis for spin-S" dimer antiferromagnets 
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The mean-field triplon analysis is developed for spin-S quantum antiferromagnets with dimerized 
ground states. For the spin-1/2 case, it reduces to the well known bond-operator mean-field theory. 
It is applied to a coupled-dimer model on square lattice, and to a model on honeycomb lattice with 
spontaneous dimerization in the ground state. Different phases in the ground state are investigated 
as a function of spin. It is found that under suitable conditions (such as strong frustration) a 
quantum ground state (dimerized singlet phase in the present study) can survive even in the limit 
5* — ^ oo. Two quick extensions of this representation are also presented. In one case, it is extended to 
include the quintet states. In another, a similar representation is worked out on a square plaquette. 
A convenient procedure for evaluating the total-spin eigenstates for a pair of quantum spins is 
presented in the appendix. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.50.Ee, 75.40.Mg 



I. INTRODUCTION 

Antiferromagnetically interacting spins arc sensitive to 
both quantum mechanics and frustration^'^. Therefore, 
the quantum antiferromagnets with frustration can re- 
alize interesting non-magnetic ground states such as the 
dimerized singlet (valence bond) states"^"^ , plaquette sin- 
glets^''', spin-liquids®"^° etc., apart from having an anti- 
ferromagnetically ordered ground state (say, Neel type). 
In the present study, we are concerned with those sys- 
tems in which the ground state has a dimer order, either 
spontaneous or given. There is an ever-increasing num- 
ber of quantum antiferromagnetic (AF) materials which 
exhibit, or seem to exhibit, dimerization physics at low 
temperatures^'-^^"-"^^. Such systems are typically charac- 
terized by an energy gap to spin excitations, thereby 
showing sharp drop in the magnetic susceptibility as the 
temperature is lowered below a certain temperature char- 
acteristic of the interaction between spins. 

We are presently interested in generic theoretical ques- 
tions concerning the instability of a dimer singlet ground 
state towards AF ordering as competing interactions in 
a system are varied. For a spin-1/2 system, such investi- 
gations at the simplest level can be conveniently carried 
out by doing a mean-field triplon analysis with respect 
to a dimer ground state. A triplon is a triplet excita- 
tion residing on a bond (dimer), and dispersing from one 
bond to another under the exchange interactions present 
in the system. An energy gap in the triplon dispersion 
implies a stable dimer phase, while gaplessness signifies 
an AF order of some kind (that depends upon the disper- 
sion). The underlying formulation is facilitated by what 
is called the bond-operator representation of the spin op- 
erators^^'^^. This approach has been successfully applied 
to many different spin-1/2 systems^^"^°. Subsequent to 
this, the bond-operator method has also been developed 
for spin-1 dimer problems^ Moreover, the triplon 
analysis has also been suitably extended to the square- 
plaquette problems (spin-1/2 case)^^. However, no such 
formulation exists for a general spin-S* dimer problem. 



In the present work, we precisely set out to achieve this 
objective. That is, to derive the bond-operator represen- 
tation for spin-S" operators (in Sec. II), and to do the 
mean-field triplon analysis for some model systems of in- 
terest (in Sees. Ill and IV). 

Obviously, it is impractical to be working with all the 
(25 -I- 1)^ states of a spin-S dimer. We therefore restrict 
our analysis to the subspace of singlet and triplet states 
only. Apart from simplifying our labor, which it does, it 
is enough for a primary discussion of the problem. Here, 
we adopt a simple working philosophy that, for an antifer- 
romagnetic spin-S dimer problem, the triplet excitations 
are the principal cause of instability (if it occurs) of a sin- 
glet phase, as the higher spin excitations are further up 
in energy and hence irrelevant for an effective low energy 
description. It is implicit in our discussion that a given 
system only has exchange interactions. The treatment, 
however, will have to be extended to include quintet or 
higher total-spin states, if the single-ion anisotropy ef- 
fects are present. While it is a real concern, presently we 
focus only on developing the triplon analysis for dimer- 
ized spin-S* quantum antiferromagnets. As an interesting 
byproduct of this exercise, we have also developed a nice 
and simple method for evaluating the total-spin eigen- 
states for a pair of spin-S (angular momentum addition) 
using Schwinger boson representation. Our method re- 
sembles that of Schwinger's, but it is different in actual 
details of the procedure that generates the compound 
spin eigenstates. ^^'^^ For the benefit of readers, it is pre- 
sented in detail in the appendix. 

Of the two spin-S quantum antiferromagnets, that we 
apply this mean-field triplon analysis to, the first one in 
Sec. Ill is a coupled columnar-dimers model on square 
lattice. Different variations of this model for the spin- 
1/2 (and spin-1) case have been studied extensively for 
investigating the quantum phase transition from dimer to 
AF ordered phases^®"'^'* . Here, we investigate this tran- 
sition as a function of S, the size of spin. Interestingly, 
we find that the dimer singlet phase (a quantum me- 
chanical phase) survives even in the so-called classical 
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limit {S — ^ oo), under suitable conditions (such as strong 
frustration). For example, we find the columnar-dimer 
phase to be stable and present in the limit 5 — ^ od in a 
small range of coupling around J2/J1 = 0.5 for the Ji- 
J2 model. While other kinds of gapped phases (such as 
the Haldane phase) are likely to arise in different regions 
of the phase diagram (especially for S > 1)'^°, we can 
not study them here within a bond-operator mean-field 
theory. However, the conclusion drawn from the present 
calculation about the 5 — >■ cx) limit, which is certainly 
valid in the vicinity of strongly dimerized limit, suggests 
a generic possibility of this kind, and asks for a rethinking 
of the classical limit in quantum antiferromagnets. 

The second one is a model on honeycomb lattice, re- 
cently constructed and shown by the present author to 
have an exact triply degenerate dimer ground state for a 
certain value of the interaction parameter, for any spin. 
In Rcf. 20; we have already presented the mean-field 
triplon analysis results for the spin-1/2 case of this model. 
However, we did not know then, how to do it for S > 1/2. 
This, in fact, was our original motivation behind develop- 
ing the triplon analysis for spin-S* dimer problems. From 
the mean-field triplon analysis of this model in Sec. IV, 
for 5 > 1, we find the dimer phase giving way to the 
Neel ordered AF phase as soon as one moves away from 
the point of exact dimer ground state. It is unlike the 
spin-1/2 case where the dimer phase was found to survive 
over a finite range of coupling. 

In Sec. V, we present two straightforward extensions 
of the spin-S bond-operator representation. First, we ex- 
tend it to include the quintet states. Next, we work out 
a similar representation on a square plaquette, in the re- 
stricted space of a plaquette singlet and certain low-lying 
triplets. Finally, we conclude with a summary. 



states. These operators are called bond-operators. 

\s):^s^\0) (la) 
IV):=iLjO) (lb) 
\Qm,)-^€jO) (Ic) 

K^,) ■■= h],„.,\0) (Id) 

Here, |0) denotes the vacuum of the bosonic Fock space. 
The completeness of the bond eigenstates implies the fol- 
lowing physical constraint on the bond-operators. 

Here, the repeated indices are summed over. 

In terms of the bond-operators, the exchange Hamil- 
tonian on the bond can be written as: 
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Furthermore, to describe the interaction between the 
spins of different bonds in the bosonic Fock space, we 
must know the spin operators in terms of the bond- 
operators. Below we develop the bond-operator repre- 
sentation for spins, which is a generalization of the well- 
known bond-operator representation for spin-1/2 opera- 
tors to the case of arbitrary spin-S*. 

In order to construct the bond-operator representation 
for the spins, we first find out the explicit forms of all 
the eigenstates on a bond. In the appendix to this pa- 
per, we have worked out an elegant procedure to write 
down the total-spin eigenstates for a pair of arbitrary 
spins. Following this approach, we can write the singlet 
wavefunction on a bond as: 



II. BOND-OPERATOR REPRESENTATION 
FOR SPIN-5 OPERATORS 



Consider the Heisenberg exchange interaction, 81-82, 
on a bond. The eigenstates, ?Tij)}, of this problem are 
such that, 8i-82|j,m,) = [-S{S + l) + \]{3 + l)]\j,m^), 
where j = 0, 1, . . . , 25 is the total-spin quantum number 
of the two spins. For a given j, the eigenvalues mj of 
operator (81 -I- 82)2 are given by m,- = — j, — j + 1, . . . , j. 
Therefore, the bond eigenstate for a given j is {2j + 1)- 
fold degenerate. It is a singlet for j = 0, triplet for j = 1, 
quintet for j ~ 2, and so on. We denote the singlet state 
as |s), and the triplets as \t^J where to^ = 0,±1. The 
quintet states are denoted as \q^^) where = 0, ±1, ±2. 
The eigenstates for j > 2 may in general be denoted as 
\hj^mj)- Below we define the bosonic creation operators, 
s^,il^^, qj^^ and h] corresponding to the bond eigen- 
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Here, the state, \S—m, —S+m), denotes a product-state, 
IS", S — m) ® \S, —S + m), of the two spins of a bond. For 
the derivation of Eq. (4), refer to Proposition 1 in the 
appendix, and also see Fig. 5. We can similarly write the 
three triplet states as: 
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where 1 in the above equation denotes = — 1 (we wiU 
sometime denote negative integers as integers with a bar), 
and the normahzation, Nt = 2S{S + 1){2S + l)/3. Refer 
to Eqs. (A. 9a) to (A. 9c) and Fig. 5 for the derivation of 
the triplet states. Below we also write the quintet states 
for 1712 = 0: 1 ^-iid 2 (refer to Fig. 6 and related discussion 
for details). 
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(6c) 
The 



2S{S + 1){2S - 1)(25 + 1)(25 + 3)/15. 
state \q^) can be generated from \q-^) by changing \S — 
m, —S + m + 1) to |5 — m — 1, —S + m) (that is, Siz O 
— S2z)- Similarly, for 192)- We can also evaluate the states 
for higher j values following the Propositions 5 and 5* in 
the appendix. 

As emphasized earlier, the present discussion will be 
restricted to the subspace of singlet and triplet states 
only. By computing the matrix elements of Si and S2 
in this restricted subspace, we derive the following bond- 
operator representation for spin-S* operators. 



S{S + 1) 
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2 \- -a ' -a~/ 2^aPyipi^ (7b) 

Here, a = x,y, z (for three components of the spin oper- 
ators), and the same for j3 and 7. The ea/B-y denotes the 
totally antisymmetric tensor. Moreover, t J. — --i=(t| — tj), 

Py ~ -^({t + f|) and il ~ fj. Since it is convenient 

to write the bond-operator representation, Eqs. (7a) and 
(7b), in terms of P^, we also write the constraint, Eq. (2), 
and the bond Hamiltonian, Eq. (3), using the same, re- 
placing <J„^C^ by PJ^. 

While our representation is valid for arbitrary S, it 
is obviously approximate for 5 > 1, because it is con- 
structed in a restricted subspace, ignoring the contribu- 
tions from j > 2 states. However, it is exact for S — 1/2, 
and correctly reproduces the known representation^^. As 
briefly discussed in the introduction, for doing a simple 
stability analysis of the dimer phase of a spin-S* quantum 
antiferromagnet, it would suffice to know the dynamics 
of triplet excitations, except when it may be necessary to 
consider higher spin states. 

In the following sections, we do mean-field triplon 
analysis of two different spin-S* models using this bond- 
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FIG. 1. Quantum spin-S coupled dimer model. As in Eq. 8, 
the exchange interactions are; thick bonds= J, thin horizontal 
lines=77a;J, not-so-thin vertical lines^r^^J, and thin dashed 
lines=77'J, where J, rjx, % and 77' > 0. Also shown are the 
primitive translations, and the spin labehng on a dimer. 



operator representation. The first one is a model of cou- 
pled columnar-dimers on square lattice. It reduces to 
many different models of interest such as J1-J2 model, 
the coupled ladders and so on. The second model is de- 
fined on the honeycomb lattice. It admits an exact dimer 
ground state for arbitrary S, and is expected to undergo 
a transition to the Neel ordered phase away from the 
exactly solvable case. 



III. COUPLED DIMERS ON SQUARE LATTICE 

Consider a spin-S* quantum antiferromagnet of inter- 
acting dimers on square lattice as shown in Fig. 1. The 
arrangement of dimers is taken to be columnar because it 
occurs in the disordered ground state of the Ji- J2 model, 
which is a special case of the present model. However, 
one may also consider other arrangements, if one wants. 
The Hamiltonian of this model is written below. 



+rjy (Sr^l ■ Sr+aji,l + Sr,2 ' ^r+ay,2) 
±ay 

+V' (Sr,l • Sr+i,2 + Sr,2 ' Sr+2ax+S,l] 



(8) 



Here, 1 and 2 denote the two spins of a dimer; r denotes 
the position of a dimer, and is summed over all dimers; 
J is the intra-dimer antifcrromagnctic exchange; various 
inter-dimer spin-exchange interactions (relative to J) are 
given by < r]x,r]y,r]' < 1. Refer to Fig. 1 for clarifi- 
cation. Physically, the single-ion anisotropy is also ex- 
pected to be present for 5 > 1 (except when the ground 
state of the ion is an orbital singlet). However, presently 
we work without such anisotropy effects. Our objective 
here is to apply the bond-operator representation of the 
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previous section to study the stability of a dimer phase in 
a simple model quantum antiferromagnet, and not study 
any particular physical system. The model given by Hi 
reduces to the following simpler models of interest for 
certain special choices of the interactions. 

1. ?/ = 0. It presents an unfrustrated model which 
interpolates between coupled dimer chains and 
coupled ladders, passing through a square lattice 
model at rj^ = rjy. For r^^^ ~ 0, it describes a set 
of decoupled (or weakly coupled) two-leg ladders, 
which is physically the case of a family of ladder 
compounds Na2T2(C204)3(H20)2, where T=Ni, 
Co, Fe and Mn in the increasing order of spin from 
5 = 1 to 5/2. As noted earlier, one must also 
consider the single-ion anisotropy in real materials, 
which has been ignored presently. A realistic calcu- 
lation for this family of ladder compounds will be 
discussed elsewhere. For rjy = 0, it describes a set 
of decoupled dimerized spin chains. 

2. ry' 7^ and Tjx = rjy = rj. It is a frustrated model 
of coupled columnar dimers on square lattice which 
reduces to the well-known J1-J2 model when rj ~ 1. 

Below the two cases are studied within a triplon mean- 
field theory of Hi. This is however not an exhaustive 
discussion of the problem, as it does not address other 
kinds of gapped phases that may arise for larger values 
of 5*. In any case, let us see what we learn from this 
simple stability analysis of the dimer phase. 



A. Mean-field triplon analysis 

Imagine a special limit of Hi in which all couplings, 
except the intra-dimer exchange J, were zero. Then, the 
spins on each thick bond (dimer) in Fig. 1 would form an 
exact singlet in the ground state. Moreover, the elemen- 
tary excitation in this case would correspond to creating a 
triplet on it. Since a triplet state on a bond costs an extra 
energy J, there is an energy gap to such excitations. Be- 
sides, these triplet excitations are localized because of the 
absence of intcr-dimcr interactions in this special case of 
independent dimers. This limit presents an idealized ver- 
sion of what is called a spin- gapped dimer phase in quan- 
tum antiferromagnets. In general, the inter-dimer spin 
interactions are non-zero, and the triplet excitations dis- 
perse, thereby lowering the spin-gap. A dispersing bond- 
triplet is often called a triplon. As long as the triplon 
spin-gap is non-zero, the dimer phase is stable against 
these excitations, and such a ground state will have zero 
magnetic moment. For some values of the competing in- 
teractions in a problem, the spin-gap may however close. 
This marks the onset of a quantum phase transition from 
the gapped dimer phase to a gapless ordered AF phase. 
We study such quantum phase transitions in Hi within 
a simple mean-field theory using bond operators. 

The key steps of a mean-field triplon analysis are the 
following. First, identify a configuration of the singlet 



forming dimers as expected in the ground state. In the 
present case, by construction, the preferred dimers are 
the tick bonds in Fig. 1. Using the bond-operator rep- 
resentation for the spins on each dimer, rewrite the spin 
Hamiltonian in terms of the bond operators, including 
the constraint by means of a Lagrange multiplier. Now 
replace the singlet bond-operators, s and s\ on every 
dimer by a mean-field, s. This results in a model of inter- 
acting triplons (on a mean-field singlet background given 
by s). To make the problem tractable, ignore the trilpon- 
triplon interaction (similar to the spin- wave analysis). 
The last two steps essentially amount to writing the spins 
on a dimer as: Sia = -820 ~ syJS{S + + i^), 

and replacing s^s by in Eqs. (2) and (3). As an ad- 
ditional simplification, we satisfy the bond-operator con- 
straint only globally. This gives us a bilinear Hamiltonian 
of triplons which can be studied fairly straightforwardly. 
In the present formulation, the quintets and higher bond 
eigncstatcs for S* > 1 arc dispersion-less higher energy ex- 
citations, and will play no role in determining the ground 
state properties. 

Applying the above prescription to Hi gives the fol- 
lowing mean-field triplon Hamiltonian. 



Hi. 



mf 



J-JSiS 



2 ^ 



[A - s^S{S + l)a] ( 
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j=2S 

EEE 

k j>2 m.j 



^k,j,mj^kj,mj (9) 



Here, Nd is the number of dimers, and = 2Jek/3. 
Moreover, = r]^ cos {2kxa) + 2(77' — r]y) cos {kyu) + 
277' cos {2kxa) cos [kyo). 

The triplon part of the mean- field Hamiltonian, H-^ 
can be diagonalized using Bogoliubov transformation. 
Define the triplon quasi-particle operators, 7i(.qj such that 



cosh6'ka7i(;Q + sinh 6'kc<7lkQj ^i^d demand that the 
triplon terms in Eq. (9) be diagonal in 7kcf This is 
achieved for tanh20ka = s2S'(S' + l)^k/[A-s25'(5'+l)^k], 
giving the following diagonal mean-field Hamiltonian. 
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(10) 



Here, ep = J - JSjS + 1) - (5 A/2) + s'^{X - J), and 
i?k = ^A[A — 2s^S{S + l)Ck] is the triplon dispersion. 
The ground state in this mean-field theory is given by 
the vacuum of the triplon quasi-particlcs, and of the ex- 
citations for j > 2. The ground state energy per dimer, 
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Eg, of the i^i is given by 



3 ^ 

eg[X,s^] = eo + ^Xl^k- 



(11) 



Minimizing with respect to A and gives the fol- 
lowing mean-field equations. 



5 _3_ V- A 
k 



3AS'(S' + 1) 

k 



(12a) 
(12b) 



The physical solution corresponds to solving these equa- 
tions self-consistently for and A. This we will do sepa- 
rately for different cases of the model. One can calculate 
the spin-gap, and also the magnetic moment in the or- 
dered phase, by solving these equations in the entire pa- 
rameter space. The staggered magnetic moment in the 
ordered AF phase is given by = 2sy/S{S + l)nc/3, 
where ric is the triplon condensate density in the gapless 
phase^^. Presently, we only compute the phase bound- 
aries between the dimer and the magnetically ordered 
phases. This is done by tracing the closing of the triplon 
gap. The wavevector Q, at which the bottom of the dis- 
persion touches zero {Eq = 0), decides the magnetic or- 
der in the AF phase. The vanishing triplon gap also fixes 
A as A* = 2s'^S{S + 1)^q. After a few steps of algebra on 
Eqs. (12a) and (12b), we get the following equation for 
the phase boundary between the columnar dimer phase 
and the Q-ordered antiferromagnetic phase in the space 
of coupling parameters. 



5- 



E 



eq - ek 



1.5 



S{S + 1) 



(13) 



In the above equation, the spin appears as S{S + 1) only 
on the right-hand-side of the equality and all the cou- 
plings are in the expression on the left-hand-side. It im- 
plies that the phase boundaries for different spins will 
collapse onto a single boundary surface in the space of 
couplings rescaled by a factor of S{S + 1) [for example, 
r]xS{S + 1) as so on]. Moreover, we can access the so- 
called "classical" limit by making the right-hand-side of 
the equality in Eq. (13) zero (that is, S — >■ oo). Below we 
will see that even in the classical limit, one finds a finite 
region of phase diagram in which the quantum mechani- 
cal singlet dimer phase survives! This seems to happens 
when the frustration is high, or when the problem is suf- 
ficiently quasi one dimensional. 



B. Calculations and discussion 

In all our calculations, J = 1 sets the unit of energy. 
Below we discuss two special cases of Hj . In the first case, 
we set 7]' = 0. This is a model of coupled two-leg ladders 
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FIG. 2. The mean-field quantum phase diagram of Hj for 
case-1 (r]' = 0). The top panel shows the phase boundaries 
between the columnar-dimer and Neel ordered phases for dif- 
ferent spins. For a given S, the region between the axes and 
the phase boundary is the spin-gapped dimer phase, and on 
the other side of the boundary is the Neel ordered phase. 
Note the quantum phase boundary in the "classical" limit 
[S — oo). Surprisingly, the gapped dimer phase survives even 
in the classical limit for sufficiently weak r]x (coupled two- 
leg ladders) or r]y (coupled chains). In the plane of rescaled 
couplings {bottom panel), the quantum phase boundaries for 
different spins collapse onto a single line. 



and coupled dimerized chains interpolating between one 
another. The second case is for 77' 7^ 0, but r]x = '/y = V- 
For 77 = 1, it is the J1-J2 model. 



1. »7' = 

We find Q = (0,7r/a) in this case. Thus, it is a case 
of quantum phase transition from the columnar dimer to 
Neel ordered phase. We compute the phase boundaries 
in the rjx-rjy plane for different values of S by solving 
Eq. (13). Here, eq — rjx + 27/y. The mean-field quan- 
tum phase diagram is presented in Fig. 2. As discussed 
earlier, the phase boundaries for different spins collapse 
onto a single line in the plane of rescaled parameters, 
r]x S{S-\-l) and r]y S{S-\-l). Fig. 2 should in principle be 
complemented with other calculations for a correct pic- 
ture in the strongly anisotropic weak dimerization cases 
(for example, to have a Haldane phase for higher spins, 
and the like) . But presently, we discuss a few interesting 
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things about this mean-field phase diagram. 

First about the model corresponding to r]x = Vy — Vy ^ 
popular fully two-dimensional case on a dimerized square 
lattice. For a finite value of 5', there occurs a quantum 
phase transition from the dimer to Neel ordered phase at 
some non-zero value of 77 = 7^*. Wc find that rj* = s(s+'i) ' 
In the limit 5 — > 00, ry* goes to zero however. That is, in 
the classical limit of this case, the ground state is Neel 
ordered even for an infinitesimally small inter-dimer cou- 
pling rj. This conforms to the usual expectations in the 
classical limit. Moreover, the spin dependence of critical 
77* ^ s{s+i) ^ simple analytical confirmation of a sug- 
gestion from the numerical studies of a similar model. '^'^ 
Besides the qualitative agreement, the mean-field calcu- 
lation overestimates the dimer phase in the present case. 
For the spin-1/2 case, 77* = 0.523 from quantum Monte 
Carlo simulations'^^''^^, 0.535 from spin- wave analysis'^^ 
and 0.54 from dimer series expansion^® as compared to 
0.62 from the present calculation. 

More importantly, we want to take note of the behav- 
ior in the extremely large S limit of the general case (that 
is, r]x 7^ Vy)- Look at the phase boundary for S = 00, in 
the top panel of Fig. 2, given by the lines: ijy = 0.022277a; 
and rjy = 11. 288772:. Here, we find two disjoint regions 
(one bounded by the lines: r]y = and r]y = 0.022277^;, 
and the other by rjy = 11.288?7j; and 77^. = 0) of the dimer 
phase (which is a quantum mechanical state with zero 
magnetic moment) existing even in this classical limit. 
This is a striking deviation from the normally expected 
behavior in the limit S 00. The two regions can be 
viewed as corresponding to the quasi- Id cases of the cou- 
pled dimerized chains and the coupled two-leg ladders, 
respectively. It seems that the strong spatial anisotropy 
in the (dimerized) lattice helps the dimer state to survive 
even when the spins are very large. Below we will see the 
same behavior also recurring in the highly frustrated sit- 
uations of a fully two-dimensional case (the present case 
of 77' = is not frustrated, but r/' 7^ in the following 
subsection is). While a case like 77^; ^ 1 and 77^ ^ 1 is 
known to be more sensitive towards Neel ordering than 
what the present calculation suggests 29,30,35^ believe 
the existence of a dimer phase for — > 00 is very likely 
to come true in more accurate numerical calculations for 
the strongly dimerized cases (77^;, 77^ ^ 1). 



2. 77' / and ^ Vy = V 

As noted above, this is the case of a frustrated two- 
dimensional model on a dimerized square lattice. In this 
case, there are two different choices of Q. For weaker 
77', Q = (0,7r/a), and for stronger 77', it is Q = (0,0). 
While the former corresponds to having Neel order in the 
ground state, the latter gives collinear order (in which 
the magnetic moments are parallel, for the spins along 
the 7/-direction in Fig. 1, and anti-parallel along the x- 
direction). From Eq. (13), we calculate the quantum 
phase diagram shown in Fig. 3. Here, eg = 377 — 477' 
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FIG. 3. The mean-field quantum phase diagram of Hj for 
case-2 (77' 7^ and rj^ = r]y = rf). Here, 77 = 1 corresponds to 
the J1-J2 model. Top panel. For a given S, the region above 
the corresponding upper phase boundary, and bounded by the 
axes, is the coUinear phase, while that below the lower phase 
boundary is the Neel phase. In between the two transition 
lines lies the columnar dimer phase. For S — 00, the phase 
boundaries are given by the equations, 77' — 0.40877 and rj' = 
0.52577. Bottom panel. The phase boundaries of for different 
spins collapse onto two lines for two different transitions. 



for Q = (0, n/a), and 477' - 77 for Q = (0, 0). For 77' 0, 
the quantum critical point for different spins is given by 
77* = s{s+i) i^^^*^ previous case), and it is 

'/* ^ c.?oX^i for 77 = 0. The phase boundaries for differ- 



S{S+1) 

ent spins collapse to a single line for the dimer to Neel 
transition, and similarly for the dimer to collinear tran- 
sition. 

The dimer phase survives again in the limit S* — > 00 
when the frustration is strong. To discuss this point, 
consider 77 = 1 case. It corresponds to the J1-J2 model. 
In the present notation, Ji = J and J2 = rj' J . In the 
classical version of this model, 77' = 0.5 is the transi- 
tion point between the Nccl and collinear ordered ground 
states'^^. It is also the point of infinitely degenerate clas- 
sical ground state manifold, and hence of very high frus- 
tration. In the quantum case, for spin-1/2 specifically, it 
is known from many numerical studies that there exists 
a quantum disordered spin-gapped state (most likely a 
columnar dimer state) in a small range of -4 < 77' < .6 
around the 0.5 point^'^"'"'. This is about .19 < 77' < .61 
from the triplon mean-field calculation^^'^^. Below and 
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above this range, one finds the Neel and colhnear ordered 
ground state, respectively, as in the classical case. 

Interestingly, even when S is arbitrarily large, wc find 
the dimer phase to be stable in a small window of 77' 
around 0.5. For 5 = 00 , this range is 0.41 ^ rj' < 0.53 
at = 1. Away from 77 = 1, the region of dimer phase 
is bounded by the lines, 77' = 0.408r/ and j/ = 0.525ry. 
Furthermore, it shrinks smoothly as one moves towards 
?7 = 0. The "quantum region" of the phase diagram man- 
ages to survive in the classical limit seemingly because 
of the strong frustration. While there may be concerns 
about the bond-operator mean-field theory overestimat- 
ing the dimer region (as it does for spin-1/2 case), a suf- 
ficient amount of frustration may always help a quantum 
state. Hence, we have a reasonable qualitative finding 
which needs to investigated further. Besides, it should 
be asked afresh, " is 5 — > 00 necessarily classical?" . 

IV. A MODEL ON HONEYCOMB LATTICE 




FIG. 4. The model of Eq. 14. The lines indicate the nearest 
neighbor Heisenberg exchange, J. The multiple spin-exchange 
proportional to K is represented by a hexagon itself, with six 
spins labeled as 1 to 6. This multiple spin-exchange inter- 
action is present on every hexagonal plaquette of the honey- 
comb. Moreover, J, K > 0. 



We now investigate a quantum spin-S* model on hon- 
eycomb lattice given by the following Hamiltonian^*^ . 



Hjj rJ ^ ^ Sj-Sj./ 

(r,r'> 



K 



Elc'2 02 q2 I q2 c;2 q2 1 
L^12'->34'->56 + '='23''45''6lJ 





(14) 



Here, J is the nearest neighbor Heisenberg interaction, 
and K denotes the strength of a multiple-spin-exchange 
interaction generated by the product of pairwise total- 
spins of three pairs of neighboring spins on a hexagonal 
plaquette. The six spins on a hexagonal plaquette are 
labeled as 1 to 6 (see Fig. 4). In the second term of Hn, 
Sfj = (Si + Sj)^. The interaction parameters J and K 
are taken to be positive. 

An important feature of this model is that it has an 
exact triply degenerate dimerized singlet ground state for 
J = for any value of 5. It presents an example of spon- 
taneous dimerization in the ground state as Hu itself has 
no preferred dimer order (unlike Hi in the previous sec- 
tion) . One of these dimer states corresponds to forming 
a singlet on every vertical nearest-neighbor bond of the 
honeycomb lattice in Fig. 4. The other two states are 
generated from the first one by making ±27r/3 rotation 
of the lattice. For J > 0, and for spin-1/2, we had earlier 
performed a triplon mean-field calculation to investigate 
the transition from dimer to Neel order in the ground 
state^*^. A similar analysis was desired for higher spins, 
but could not be done at that point due to the lack of 
spin-5 bond-operator representation. It was our original 
motivation for developing the bond-operator mean-field 
theory for arbitrary spins. Having achieved this objective 
in Sec. II, we can now do a triplon mean-field theory for 
iJii, exactly in the same way as done for Hj, by taking 
the configuration of vertical dimcrs on honeycomb lattice 
as a reference state. 



The mean-field triplon Hamiltonian in this case also 
looks the same as in Eq. (9). The diagnoalized mean- 
field Hamiltonian for Hu can therefore be written as: 

1 



-|-(localized higher j dimer states). 



(15) 



IX+s^{X~J~2K[SiS + l)]^) + J[l~S{S + 



where eo ^ 

1)] + 2K[S{S+1)]^, a.nd^i,^lS{S + l){J + 2K[S{S + 
1)]2(1— .s2)}efc. Here, = cos 27rfc2+cos 27r(fci — ^2)7 and 
k is defined as k = 27r(fciGi 4- ^202) where Gi and G2 
are the primitive reciprocal lattice vectors, and fci,/e2 G 
[0, 1]. The self-consistent equations for this problem are 
given by 



3 
2N'd 



E 



A 



VA(A - 2s2^k) 



and 

A = J + 2K[S{S+l)]'^ 



3 

2iVd 



E 



2 a«k 



v/A(A - 2s-2ek) 



(16) 



(17) 



The closing of the triplon gap, which marks the insta- 
bility of the dimer phase to a magnetically ordered phase, 
fixes A as A* = 2s^£^q, where Q = (0,0) in the present 
case. At this quantum critical point, we get 



■E 



2 2iV.^yT-f 



0.817, and 



(18a) 



2[S{S+1)]'' 



2s' 



3* 



S(S+1) 



8 h2 
3* 



0.71 



1.524- 



1 

S(S+1) 



a (18b) 
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where u 



N, 



0.655. Equation (18b) is a 



closed form expression for the critical J/K, denoted as 
as a function of S. For the dimer order to become 
unstable to Ncel order in this mean-field theory, J/K 
must be greater than C,*g for a given S. We find Q'^^^ ~ 
3.067, which agrees with our earlier calculation for the 
spin-1/2 case of this model. Next we find Ci — —2.454, 
(■3^2 — —12.247, and so on. For positive J and K, we 
therefore conclude that the mean-field triplon calculation 
predicts a Neel ordered ground state for Hu for 5 > 1 
for any non-zero value of J / K . 

We know for sure that the dimer ground state is exact 
for J = 0. For the spin-1/2 case, the triplon analysis 
predicts that the dimer state will give way to the Neel 
state only when J/K is sufficiently strong. However, 
for spin-1 and higher, it seems to happen for arbitrar- 
ily small J. At this point, it is important to note the 
following. While deriving the mean-field triplon Hamil- 
tonian for Hn, following the steps outlined in the pre- 
vious section, we end up having no contributions from 
the six-spin terms of the form (Si • S2)(S3 • S4)(S5 • Sg). 
For iS* = 1/2, it does not seem to affect the dimer phase 
for small J as the mean-field theory suggests. We have 
some evidence of this from a numerical calculation in our 
earlier work on this model^°. However, it is not clear 
as to how the absence of contribution from the six-spin 
terms in the present mean-field calculation will affect the 
case of spin-1 and higher. May be, in a renormalized 
triplon analysis of Hu, one gets the dimer phase over 
a small but finite range of J/K ior S > 1. One such 
calculation is done by writing (Si • S2)(S3 • S4)(S5 • Sg) 
as -[Sis + l)]^s^x [S3 • S4 + S5 • Se - 2x^(5 + 1)], and 
similarly for (S2 • S3)(S4 • S5)(S6 - Si). Here, the expecta- 
tion values (S3 • S4), (S5 • Se), (Sg • Si) and (S2 • S3) are 
all taken to be equal to S{S + l)x- It gives the following 
critical value of J/K. 



|^2[5(5- 



0.856 



1.524^ 



- 1 



(19) 



S{S+1) 

It is similar to Eq. (18b), except the numerator inside 
the square-brackets in different, which only slightly in- 
creases the value of (g. But the qualitative conclusion 
remains the same. That is, the Neel order sets in for ar- 
bitrarily small J for > 1. Well, this is the result from 
triplon mean-field calculation. Alternative calculations 
are needed to resolve this conclusively. 



V. EXTENSIONS OF THE REPRESENTATION 

Below we present two immediate extensions of the 
bond-operator representation derived in Sec. II. First, we 
go beyond singlet and triplets to include quintet states 
on a bond. In the second case, we derive a similar rep- 
resentation on square plaquette in terms of the plaque- 
tte bosons, which turns out to be an easy extension of 
Eqs. (7a) and (7b) to a plaquette problem. 



A. Including quintets on a bond 

The bond-operator representation derived in Sec. II is 
in the subspace of singlet and triplet states only. With 
some labor, we can extend this to include the quin- 
tets, knowing how to systematically construct the bond- 
eigenstates (see Appendix). Including higher states is 
possible, but it requires even more effort, and will not be 
considered presently. The bond operator representation 
including quintets is written below. 



-[i\q, + qli,+i{q-^+qli-. 



qUi - qUi 




(20) 




1 , 



7! (^1^0 - qlk) 



(qUi + qlh) 



(21) 



Here, N,, = 25' -I- 1, = 25(5' + l)N^/3 and N, = 
Nf(25— l)(25 + 3)/5, are the normalization constants for 
the singlet, triplet and quitet states respectively. More- 
over, in the notation ±, the '+' corresponds to Si and 
'-' to S2. 

Equations (20) and (21) are exact for spin-1 case^^, 
and reduce to the representation for spin-1/2 operators 
by dropping the terms involving quintets. Note that the 
coefficients of the terms mixing singlet with triplets and 
triplets with quintets scale as 5' for large 5. While the 
strengths of different mixing terms grow similarly as S 
grows large, the hierarchy of mixing suggests that for a 
spin-gapped phase in a system of exchange-interacting 
quatum spins, the triplon analysis is a minimal reason- 
able thing to do. It is because the condensation of "quin- 
tons" is facilitated only by that of the triplons. In a 
gapped phase where triplons have not condensed, it is 
unlikely that the quintons will condense. Therefore, it 
seems okay to ignore the quintet states to first approxi- 
mation. It will not be the same however if we take into 
account the single anisotropy effects like (5f 2)'^- In this 
case, the singlet state will directly mix with quintets, 
and therefore, it will be better to work with Eqs. (20) 
and (21) instead of Eqs. (7a) and (7b). This combined 
"tiplon-quinton" analysis will be useful in investigating 
the infiuence of single-ion anisotropy on the stability of 
a dimer phase, and on its existence in the limit 5 — >■ 00. 
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B. Representation on a square plaquette 



Consider a spin- 5* problem on a single square plaquette 
given by the Hamiltonian: Hgp ~ J(Si • S2 + S2 ■ S3 + S3 • 
S4-|-S4-Si)-|- J'(Si •S3-I-S2-S4), where J is the exehange 
interaction along the edges of the square and J' is the in- 
teraction along the diagonals. The subscript sp stands for 
square plaquette. This problem can be solved by rewrit- 
ingit as: H,j, = ^S^ - ^^(Sfg + Si4) - 2J'5(5 + 1), 
where Stot = Si +S2 + S3 + S4 is the total spin of the pla- 
quette, and Si3 = S1+S3 and S24 ~ S2 + S4 are the total 
spins on the two diagonals. The eigenstates of this prob- 
lem are completely specified by three quantum numbers: 
the total spin of the plaquette, and the two diago- 
nal spins, ji3 and ^24, with eigenvalues, ji3, J24) = 
4j-(j + l)-^^bi3(ji3 + l)+j24(j24 + l)]-2J'^(5 + l). 
Given that we are interested in antiferromagnetic inter- 
actions (J, J' > 0), let us figure out the possible ground 
states, and derive a bosonic representation for spin oper- 
ators, considering only the lowest energy excitations. 

Since J > 0, for a given J13 and j24, the E^p would be 
lowest for the smallest value of j. Moreover, when J > J', 
the ground state of Hgp is given by: ji3 = j24 = 25* and 
j = 0, and for J < J' , it corresponds to ji3 = ^24 = 0. 
The latter is a case of dimer ground state, in which 
the two diagonal bonds separately become singlet, and 
j is trivially zero. The elementary excitations in this 
case would just correspond to making a diagonal bond 
a triplet. In short, for J < J', the bond-operator rep- 
resentation of Eqs. (7a) and (7b) is applicable as it is. 
However, for J > J', the ground state is a true plaquette- 
singlet, involving all four spins. Therefore, we must sep- 
arately find out a representation of the spin operators 
in terms of this plaquette-singlet and the corresponding 
plaquette-triplet excitations of the elementary kind. 

For J > J' , the ground state lies in the sector given 
by .713 = J24 — 25*. The plaquette states in this sector, 
for different values of j, are the compound eigenstates of 
two spins of size 2S. That is, in Eqs. (4) and (5a-5c), 
replace S by 25. This immediately suggests that the di- 
agonal spins, S13 and S24, are represented by Eqs. (7a) 
and (7b) with S written as 25', where s and ta opera- 
tors are now the bosons corresponding to the plaquette 
singlet, \j — 0;ji3 = 25", ^24 — 25*), and triplet states, 
\j = l;ji3 = 25*, j24 = 25), respectively. In order to find 
the representation for individual spins. Si and S3, we 
must also find Si — S3 in terms of the plaquette bosons, 
and do similarly for S2 and S4. Since we consider only 
those states given by ji3 = j24 ~ 25, the operators 
Si — S3 and S2 — S4 would be null operators in this 
restricted subspace because they change the values of ji3 
and j24- Hence, their matrix elements in the subspace 
of the ground state singlet and the lowest triplets are 
zero. It leads to the following representation of the spin 



operators on the plaquette. 
Siq ~ S3Q, 

S2a = S4Q 

~ -^M^^+lI (st£^ + £t^^) _ le„^^£^i^22b) 

The above equations are written in the standard nota- 
tion, except that the bosons are now defined on a square 
plaquette. It correctly reproduces the representation for 
spin- 1/2 case^^. While it is an approximate representa- 
tion, it provides a simple framework for discussing the 
low energy physics of a (coupled) plaquette problem for 
J sufficiently stronger than J' (and other couplings in a 
given problem). However, when J' is strong enough, the 
states from other sectors begin to compete. For example, 
on a single plaquette for J' > (1 — -^) J, the singlet state 
for J13 = ^24 = 25 — 1 becomes lower in energy than 
the triplets in the sector containing ground state. This 
renders the above representation insufficient for an effec- 
tive low energy description. It is, in any case, a useful 
representation, if considered within limits. 



VI. CONCLUSION 

To summarize, we have derived the bond-operator rep- 
resentation for spin-S dimer problems, and also worked 
out a similar representation on a square plaquette. Us- 
ing this bond-operator representation, we have done the 
mean-field triplon analysis of two model quantum anti- 
ferromagnets: 1) a coupled columnar dimers model on 
square lattice, and 2) a model on honeycomb lattice with 
spontaneous dimer order in the ground state. Through 
this mean-field calculation, we have studied the quantum 
phase transition from the dimer to AF ordered phases as 
a function of spin. A notable outcome of this analysis is 
that one finds the dimer phase, which is a quantum me- 
chanical phase, to exist even in the limit 5 — >■ 00, under 
the conditions of strong frustration (or spatial anisotropy 
with strong dimerization) . It suggests that the limit 
5 — >■ cx) is not necessarily "classical" , as there may not ex- 
ist any ground state with non-zero classical magnetic mo- 
ments for a system of (frustrated) quantum spins. Such 
quantum ground states are known to exist for arbitrar- 
ily large spins in specially constructed models, such as 
the Shastry-Sutherland model^ or the exactly solvable 
case of the model on honeycomb lattice in Sec. IV. How- 
ever, we believe this behavior of having quantum states 
in the so-called classical limit to occur more generically. 
The present observations offer an interesting view on the 
classical limit of frustrated quantum antiferromagnets, 
which further needs to be investigated carefully. 
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Appendix: Compounding a pair of spins 

Here, we describe an interesting approach, that we 
have developed, for adding two quantum spins using the 
Schwinger-boson representation for spin operators. It 
generates the closed-form expressions for the Clebsch- 
Gordan coefficients rather conveniently. 

Let d and b denote the Schwinger-boson operators, in 
terms of which, the operators of a spin can be written as: 
S+ = dt5, S- = S^a, and = (ata - StS)/2, subjected 
to the constraint: a^d + b^b = 2S. This is called the 
Schwinger-boson representation. Here, and Sz are the 
usual spin operators, and S is the spin quantum number. 

For a pair of spins Si and S2, define an antisymmetric 
pair-operator (also called the valence-bond operator)'''^: 



a\bl 



b\al, 



(A.l) 

(A.2a) 
(A.2b) 
(A.2c) 

where the subscripts 1 and 2 stand for the two spins (or, 
two sites of a bond) . These pair-operators will form the 
basis of our analysis for evaluating the total-spin eigen- 
states of a pair of spins. Besides, we will also use the 
following operator identities for systematic proofs. 



and three symmetric pair-operators: 

Bi = a\al 

Bl = albl + b\al 

B\ = b\bl 



(A.3) 



Here, A/i = a[ai + b\bj^, and Af2 = Ojdj -I- 62^2 ^''^ 
number operators of the Schwinger bosons on site 1 and 
2, respectively. By successive application of the above 
relation, we can show that 

A''A{A''y |ni,ri2) = {A'^y /(ril +712 + / + !) 



-A^A 



\ni,n2) 



(A.4) 



where |ni,n2) denotes a state with total number of ni 
Schwinger bosons on site 1, and n2 on site 2. The other 
useful relations are: 



1 



2Si • S2 = -A/'iAA, - A'' A 



(A.5) 



and 



(S1 + S2) 



-{JVi+M2){Mi+M2 + 2)-AU. (A.6) 



We will also need the following commutators. 



A,{a\y 


- i{a\y-% 


(A.7a) 


'A,ib\y 


- ^i{b\y-'a2 


(A.7b) 


A,{aly_ 


= -ii4y-'h 


(A.7c) 


A,{biy 




(A.7d) 



Now we are all set to construct the total-spin eigenstates, 
I J, TOj), where j denotes the total-spin quantum number, 
and nij is the quantum number for the z-component of 
the total spin for a given j. 



1. Case of equal spins: Si = S2 = S 

• Proposition 1. The normalized singlet eigenstate 
is given by: 



|j =0,m, =0) = 



1 



{2S)\y/2S- 



= {A^f\m 



Proof. It is clear that ni = 712 = 25' in the proposed 
state, (^^)^'^ |0, 0), where |0,0) denotes the Schwinger- 
boson vacuum in which ni = n2 — 0- 

Apply (Si + 82)^ on the proposed state. Using 
Eqs. (A.6) and (A.4), we find that 

(Si + S 



(y^t)2^|0,0) =0. 



Since the proposed state is annihilated by the total-spin 
operator, it is a singlet. That is, j = 0. Moreover, 

(Si + 82)2 (A^)"^^ |0, 0) = 0, because A^ changes the to- 
tal number of a-type and 6-type Schwinger bosons by the 
same amount (that is, one). Therefore, rrij = 0. Now, 
we fix the normalization. 

Let, NormpS'] (0, 0|(yl)2S (ylt^ 
malization constant. Clearly, Norm[0 
find that 



.2S 



0, 0), be the nor- 
: 1. Moreover, we 



Norm[2S'] = (25 + l).2S'.Norm[2S' - 1]. 

It is derived using AA^ = (2-F7Vi +7V2)(4+7Vi +7V2)/4- 
(Si + 82)^. This recursive relation for normalization im- 
plies, Norm[25] = {2S + l)[{2Sy.]'^ . Hence, the proof. • 

Next we workout a procedure for generating the eigen- 
states for arbitrary j'. Let us introduce a "generating" op- 
erator, B'^{^) = B\+£_Bq+£,^B\^, where ^ is just a parame- 
ter. It can also be written as: B^{^) = {a\+£,b\){a\+^b2). 

• Proposition 2. The total-spin quantum number of 



the state, = [A 



0,0), is J. 



Proof. Evaluate (Si + "^2) as follows. 

(Si+S2)2|j;C) 

= [25(25 + 1) -ylU] {A^f"^"' [6^(0]' |0,0) 
[Equation (A.4) implies the following] 

= j(j + i)|j;0 - A^^-^^^A [st(0]' |o,o) 

For |j;^) to be an eigenstate of the total-spin operator 
with eigenvalue j(j ~\- 1), the operator A must annihilate 
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|0,0). Below, wc show that it is true. 



[Eqs. (A. 7a) & (A. 7b) imply the following.] 

E cicie^+'^ ia\y-'^Aib\y^ ialy-'^ (biy^- 10, 0) 

Zl,/2=0 

" V ' 

E C/, C/^e'^+'- ia\y-'' ib\y^ Ajaly-'^ jblt' |0, 0)_ 
= 

Here, the quantities above the under-braces are zero (first 
one due to a simple cancellation, and second because the 
position of A is such that ai and bi annihilate |0, 0)). The 
relation, A [B'^{(,)y |0, 0) = 0, straightforwardly implies 
(Si + S2)2|i;0 + Hence, the proof. . 

Below we show that all the different rrij states for a 
given j are contained in |j; ^), and express themselves in 
powers of ^. Since all |j, rrij) states can be systematically 
derived from \ we call \ the "generating state". 
Note that |j; ^) is not a normalized state. 

• Proposition 3. The generating state, \ has the 
following series expansion in powers of ^. 



b;0= E e'-'-^b'^-) 



Here, |j, m^) denotes an unnormalized total-spin eigen- 
state (as compared to |j, rrij), which denotes the normal- 
ized version of |j, rrij)). 

Proof. To derive the series form of \j]^), expand 
[i3^(^)]"' in powers of ^. 

= i2J2 ciciie'+'HB\)=-'HBiy^-'HBiy^ 

ll=0l2=0 

Note that B\ contributes 1 to (Si + 82)^, Bq contributes 
zero, and contributes —1 to the same. Moreover, A^ 
adds nothing to it. Therefore, the eigenvalue of (81+82)2 
for the terms corresponding to ^'i+'s in the generating 
state, \j; ^), is equal to j — h — h- Make a transformation 
of the summation variables from (^1,^2) to {rrij, I) such 
that rrij = j — h — h and I ^ I2. Or conversely, h = 
J — rrij — I, and I2 — I, where rrij = j — t- —j and I = 



max(0, —nij 



2 



Here, [x] denotes the integer- 



valued part of X. Now, we can write: 



where 



E 

i=max(Or 



{rrij + iy.{j - - 2l)Ul 



X (^1)™^+' (4)^-™^-2' {Biy (A.8) 

The (unnormalized) generating state is thus written 
as: =T.7n,=j^^~"''\j,rnj), where 

|j;^)=6t(j,m,)(^t)2S-J|o^O). 

It has been argued above that (81 + S2)z\j,mj) = 

mj\j,mj). Therefore, {j,nij\j,m'-) = for rrij ^ m'y 
Furthermore, Proposition 1 implies 

E e'""^- [(Si + 82)^ - 3{] + 1)] il^) = 



which in turn implies (Si + 82)^!^, nij) — j(j + l)|j, rrij). 
Hence the proof. • 

The above mathematical result can be understood in 
the following way. Think of the singlet state with 25* va- 
lence bonds as a reference state, kind of a "valence-bond 
sea". Out of which, one can generate different \j,mj) 
states by removing j valence bonds, and inserting the 
same number of symmetrized bonds in a suitable way. 
This insertion is precisely given by the operator B^{j, rrij) 
defined in Eq. (A.8). For example, the triplet eigenstates 
(j = 1 and nij = 1, 0, —1) can be constructed as: 



\yi)^Bl{Mf' '|o,o) 



\l,0)ocBi{A^f' >,0) 



\l,-l)^Bl{A^f' '|0,0) 



(A.9a) 
(A.9b) 
(A.9c) 



This procedure is pictorially illustrated in Fig. 5. 

Similarly, the eigenstates for j = 2 and rrij — 2,1,0 
are given by (see Fig. 6): 



\2,l)^BlBi{A^)''-'\0,0) 



(A.lOa) 
(A.lOb) 



|2,0)(x {Bl 



The states for 



2B\bI 



{A^Y"^ '|0,0) (A.lOc) 



— 1 and —2 can be obtained by 
replacing B\ by in the equations for |2, 1) and |2,2), 
respectively. 
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Singlet State 



^[albl-b\al) |0) 



2S 



(VBs = Valence Bonds) 




■<o>- 



2S-1 
VBs 



|i,o) 



■<1>- 



2S-I 
VBs 



11,-1) 



FIG. 5. (Top) The singlet state for a pair of spin-S. It can be 
viewed of 25* valence bonds. (Bottom) The three 

triplet states are created by replacing one valence-bond by 
three different symmetric bonds. See the text for details. 



2S-2 
VBs 



|2,2) 
-©- 



2S-2 
VBs 



|2,1) 

R— ®- 

— (5)— 



2S-2 
VBs 



+ 2x 




|2,0) 



FIG. 6. The quintet states. See Eqs. (A. 10a) to (A.lOc). 



While Eq. (A. 8) can be nicely visualized, and is helpful 
in understanding the generating procedure, it is not the 
most convenient form of (j, nij ) . Below we present a 
more effective form of this operator for evaluating |j, rrij). 



B\j,Tnj 



mm{3,j-mj) 
i=max(o,-mj) 



X (al)™^+' {b\y-^'-' (at)-'--' {b\)' (A.ii) 

The above expression is derived by noting that B^{^) = 
{a\ + £,b\){a\ + £b\)- Having learnt the generating state 
description of the total-spin eigenstates, we now find the 



normalization constant, (j, | j, ) . 

• Proposition 4- The normalized total-spin eigen- 
state, with quantum numbers j and nij, is given by: 



I j> 



1 



1.7, "^j-; 



where |j, m^) are the corresponding unnormalized states 
in Proposition 3. 

Proof. Following the same steps as for the normaliza- 
tion of the singlet state in Proposition 1, we find that 

We further find that 

{o,o\m)y[B\oy\o,o)^[m+ey]'- 

Therefore, 

Moreover, = Em,=j C^*^'""'' |j>^-) (as 

implied by Proposition 3). Therefore, 

and the normalized total-spin eigenstate, \j,mj) is: 

1 



I J: 



J'- 



ni2s-jy.Jc',f+^-''c'i 



Hence, the proof. • 

Finally, we derive the explicit expression for \j, rUj) in 
terms of the actual spin states. It will give us all the 
Clcbsh-Gordan coefficients in a closed form for artbirary 
i and ruj for a given pair of spin-S". 

• Proposition 5. The total-spin eigenstates, |j,mj), 
can be written as: 

2S-\mj\ 
m=0 

\S — m, —S + m + iTLj) V < mj < j 

\S — m + rUj, —S + m) V — j < nij < 

where the Clebsch-Gordan coefficients, C[j, mj;m], are 
given by^^: 



C[j,mj; 



^2S+j+1^2j ^2S ^23 



mm(2S-j,m) 

E 

p=max{0.m+\mj\-j) 



( \pr~i2S-j f^j r-ij 
\ I ^m-p^m+\mj\-p 
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The states, \S — m, —S + m + rrij), above denote the 
product states, |51 = 5, toi = 5 — m) (g) = S,m2 = 
— S + m + rrij), of the two spins. 

Proof. Consider B^{j,mj) as given in Eq. (A. 11). We 
discuss the positive and negative values of rrij separately. 

For < rrij < j, 



Therefore, 

b^)=St(j,m,)(^t)2S-^|0,0) 

- E E(-)''^™,+.^^/cf-^'(«I)^"^-''+^^-^-^^ 

;=o p=o 

x(5t)0""b-'+p) (4)0'-'+p) (^)('+2S-j-p)|Q^o) 

x\S^j-p + l + mj,-S + j+p-l) 
[change of variable: I ^ j — rrij — I] 

j-nij 2S-] (-<3(-<3 r<2S~3 

'=0 P=0 \J'^'l+p'^l+p+rnj 

x\S - I - p,-S + l+p + nij) 

Define new variables m and p as: m = I +p, p = p. Now, 
we can write the normalized state, |j, mj), as: 



2S-mj 

E 



1 



''m-\-mj—p^p 



Ip2s+2j+ip2j ^ Ir^sr'^s 

y^-2j+l ^j+m, ™=0 Y'-^m'-^m+r: 

min(2S-j,m) 

yp=max(0,m+mj - j) 
X — 771, — S* + 771 + 777j) 

= C^[j, 7Tij; mJIS* — m, — 5 + m + mj) (A. 12) 



jn=0 



where the coefficients of linear combination, famously 
called the Clcbsch-Gordan coefficients, are given by 



C2S 



p2S+2j + lp2j p2S r2S 
^'2j + l ^j+m-j^m+mj^m 



mm(2S-j,m) 

E {~Y'~^m-p^taj+m-p'~^p^ ^ (A. 13) 

p=max(0,m+Tnj-j) 

Here, the superscript, +, indicates that it is for positive 
values of rrij . Also note the summation variable p written 
as p (it is allowed for dummy variables). We can carry out 
the same analysis for the negative values of ttt,-,-. However, 



we will infer the negative ruj states using an interesting 
argument described below. 

For —j < rrij < 0, write rrij ~ -~\mj\. Then, 

6t(j,m,) = 6t(j,-K|) 

l=\mj I 

[change of variable: 1-^1 + \mj\] 

Comparing the last line of the above equation with 
B^{j, \mj\) suggests that rrij — > —rrij is equivalent to the 

mapping: ai o 62 and bi -O- 0,2- Under this mapping, 
the valence-bond operator A is invariant. However, the 
spin quantum numbers exchange: 5*1 O 82- Moreover, 
Siz O — S22. Therefore, the total-spin eigenstate for a 
negative mj is given by 

\3,mj) = \ j,-\mj\) 

2S-\mj\ 

= 'C'^[j,\mj\;m]\S — m ~ \m.j\,—S + m) 

m=0 
2S-\mj\ 

= C^[j,mj;m]\S — m + rrij, —S + m) (A. 14) 

m=0 

where the coefficient [j, 777,-, ; 777], for negative 777^-, is 
equal to C"'"[j, |777j|; 777]. Equations (A. 12) and (A. 14) to- 
gether can be stated as follows. 

2S-\mj\ 

\S — 777, —S + 777 -I- 777j) V < 777j < j 
jS* — 777 -I- 777j, —S + m) V — j < 777j < 

Or even more compactly, 

2S-\mj\ 

b'-'^i) = E C[j, 777j;777]X 
m=0 

\S — 777 + min(0, 777^), ~S + m + max(0, 777j)(A.15) 

where C[j, 771^ ; 777] = C^[j, |777j|; 777]. Hence, the proof. • 
With some care, we can write C[j, j^ ; 777] in the follow- 
ing very compact form.^^ 

C[j,777,; 777] = (-)'»+l"bl-^- C^„^ ,C7+I'"^I X 



(J2S 
m+\mj I 



^2S+j+1^2j ^2S 



,F2[a,b, 1](A.16) 



where 3F2 [a, h, 1] is the generalized Hypergcometric func- 
tion, and a and h (two arrays of size 3 and 2, respectively) 
are given by: a = + |r77j|,r77 -I- \mj\ — 2S} and 
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b = {l + lmjl, l—j + m+\mj\}. (For the definition oipFq, 
look up in any book on special functions or Matlieniatica 
or Google.) 

We can now explicitly write down any total-spin eigen- 
state for a pair of spin-5. For example, the singlet state 
(denoted as |s)) can easily be written as: 



V25 



2S 

= y{-) 



\S-^ 



We can similarly evaluate the triplets and other higher 
spin states. Next, we generalize our method to the case 
of unequal spin c^uantum numbers, S'l and 82- 



2. Case of general 5i and S2 

• Proposition 2*. The generating state for the 

total-spin quantum number j, is given by: 



j + Sl-S2 



j + S2-Sl 



|0,0) 



where \Si ~ 82] < j < Si + S2. 

Proof. This is a generalization of Proposition 2 (hence, 
2*). Likewise, evaluate (Si + S2)^|j;0- We find that 

{Si + S2f\j;0 

= j{j + i)|j; - [B\oy |o, 0) 

Following the same steps as in Proposition 2, we can 
show that A^'+^^-^^+^A[B^^)y |0,0) = 0. Therefore, 

(Si + S2)2|j-;0 = .?(.? + i)U;0- 

The range of j is fixed by demanding that the powers 
of {a\ + £,b\), {a\ + ^b\) and in the state must 
be positive integers. It implies that j > |5i — 5*2! and 
j < Si + S2. Physically, the lower bound, | S'l — 52], is 
tied to the fact that a maximum of min(2S'i, 2S2) valence 
bonds can be made between two spins. The upper bound 
on j is fixed by the total number of Schwingcr bosons, 
2(51 -I- 5*2). Hence, the proof. •. 

• Proposition 3*. Series expansion of 



Here, \j, rrij) is an unnormalized eigenstate of (Si -|- S2)^. 

Proof. Consider {a\ + ^61)^+^1-^^(4 + fbly+^^-^^ 
first. It can be expanded as: 



= g e^-™^St(j,m,) 



where 

i+min(S2-Si,-mj) 
i=max(o,S2-Si-mj) 

The above equation is a generalized version of Eq. (A. 8). 
Now it's obvious that |j; = 'I2ni =j \ j, mj), where 



Sl+S2-j 



|0,0). 



The above result, together with Proposition 2*, further 
implies that (Si -fS2)^ I j, TOj) = j{j+l)\j,mj) . Moreover, 
Eq. (A. 17) implies (Si + 82)21.7, mj) = ™ib7™j) because 
A^ contributes nothing to (Si + 82)^. • 

• Proposition 4* . Normalized total-spin eigenstate: 



Si+S2-j^Si + S2-3 



C: 



■2S2 



-Ml my 



^{2Si)\{2S2)\Cl^ty^^^+'C%,^^ 
Proof. Calculate i'^ Proposition 4. We find, 

(j; ^b; = cl^ti'^'^' [{Si + S2 - j)f x 
{3 + Si-S2)\ {] + S2-Si)\ [i + et'- 

Moreover, {j]i\i;Cl = Em^-=i 0'. "^j li. "^j)- de- 

duced from Proposition 3*. Therefore, 

[(Si + S2- 3)f U + Si - S2)! U + S2 - Si)! 

_ {2Si)i{2S2)ici^yy^+=+'c]i^^ 

^S^^S2-j^Sl+S2-j 

Hence, the proof. • 

• Proposition 5* . The Clebsch-Gordan coefficients. 

|j, my = E^ C[j, mj\ m] |Si — m, —Si -f- m 4- my 

m—TClmin 

where 

mmi-n = -min(0,mj -l- S2 - Si), 
-mraax = Si + S2 - max(mj, S2 - Si), 

and the Clebsch-Gordan coefficients, C[j, mj;m], are 
given by 



■2Si 



2S2 



Si+S2-3^Sx+S2-j 



^Si + S2+i+1^2j ^2S2 -^2Si 
'-^2j + l ^ 3^mj^m+mjJrS2-S\^ 

\p/^i-i-Si-S2^i+S2-Si ^Si+S2-j 



15 



where 

Pmin ~ inax[0, m + niax(TOj, 5*2 — Si) — j], and 
Pmax = min[S'i + S2 - j, rn + min(0, mj + 82-81)]. 

Proof. One can get it from Propostions 3* and 4*, 



by carefully doing a few steps of algebra, similar to that 
in Proposition 5. Note that this proposition correctly 
reproduces Proposition 5 for 5*1 = S'2 = S.» 

This completes our description of a Scliwinger-boson 
based method of constructing the total-spin eigenstates 
for a pair of quantum spins. 
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